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Abstract 

We investigate mesons in the bosonized massive Schwinger model in the 

hght-front Tamm-Dancoff approximation in the strong couphng region. We 

confirm that the three-meson bound state has a few percent fermion six-body 

component in the strong coupling region when expressed in terms of fermion 

variables, consistent with our previous calculations. We also discuss some 

qualitative features of the three-meson bound state based on the information 

about the wave function. 
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I. INTRODUCTION 



Recently there has been increasing interest in hght-front field theory [^. In particular 
the light-front Tamm-Dancoff (LFTD) approximation p| has proven to be very powerful, as 
a nonperturbative approach to the relativistic bound state problem. It has been successfully 
applied to two-dimensional models as well as to four- dimensional Yukawa theory 0]. It 
is important to note that this new approach not only reproduced known results correctly, 
but also brought us new results, which have never been obtained by other methods 0,0. 
(See also Refs. for some of such "new" results in the discretized light-cone quantization 
(DLCQ) approach §.) 

In a previous paper we studied the massive Schwinger model in the framework of 
the LFTD approximation up to including fermion six-body states. We showed that (1) 
the two- meson bound state has a negligibly small six-body component, (2) the three-meson 
bound state does exist, and (3) the two-meson bound state is well described in terms of 
the wave function of the relative motion. Surprisingly, however, the six-body component of 
the three-meson bound state is quite small, though it is large in comparison with those of 
other states below the three-meson threshold. Typically it is at most a few percent for a 
small value of the fermion mass. Despite this small six-body component, we identified it as 
the three-meson bound state based on the following reasons: (1) Since the meson creation 
operator contains the fermion annihilation operators (see Ref. for the notation), 

A\p)= '4){k,p-k)h\k)d\p-k) 

^0 {2'K)^k{p-k) 

+ / i^{p+k,k)[h\p+k)h{k)-d\p + k)d{k)], (1.1) 

•^0 {2tx) ^k{p + k) 

the three-meson bound state (~ A^'^ |0)) naturally contains fermion two- and four-body 
components, beside the six-body component. (2) This state is charge conjugation odd, 
while a two- meson state should be charge conjugation even. 

Unfortunately, we were not able to explain why the six-body component of the three- 
meson bound state is so small. It is not the aim of this paper to do it, but to further confirm 
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our conclusion and to give some arguments by using bosonization. 

How can we further justify the identification despite the smallness of the six-body com- 
ponent of the state? As emphasized in the previous paper, we have a simple picture of the 
massive Schwinger model in the strong coupling region. Because the massless theory (the 
strong coupling limit) is a free massive boson (meson) theory, we expect that it becomes 
a weakly interacting massive boson theory once a small mass term is included (the strong 
coupling region). In addition, light-front field theory provides us with the simple vacuum. 
Actually, these two allowed us to describe the two-meson bound state in terms of the wave 
function of the relative motion §]. One may think, therefore, of constructing the wave 
function of the relative motion for the three-meson bound state in a similar way. If such a 
description is a good approximation of the state, one may justify that it is a three-meson 
state. It turns out, however, that it is almost infeasible in terms of the fermion variables 
because it is difficult to find a simple set of basis functions which satisfies all the symmetry 
properties. 

A simpler way we will follow in this paper is to consider the bosonized theory. 

Is the Tamm-Dancoff approximation good for the bosonized theory too? The meson 
has an internal structure in the sense that the wave function, ip{k,p — k) in eq. ( |1 . 1| ) , has 
nontrivial momentum dependence, which becomes negligible in the strong coupling limit, 
i.e., the meson becomes "structureless." The Fock states in the bosonized theory is that 
of this "structureless" meson. The "structureless" meson state is a good approximation in 
the strong coupling region. It means that the Tamm-Dancoff approximation is good for the 
bosonized theory in the strong coupling region. (The wave function gets the momentum 
dependence from many-body states.) In the bosonized theory, it is easy to find a set of basis 
functions which satisfies all the symmetry properties. 

In this paper, we consider the bosonized massive Schwinger model in the LFTD ap- 
proximation up to including three-boson states. We show that almost all of our previous 
calculations ||^ are consistent with the results obtained in the bosonized theory in the strong 
coupling region. In particular, we show that the three-meson bound state, which is almost 
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100% a three-boson state, has a few percent fermion six-body component. 

We are aware of the hmitation of our analysis in this paper. The hmitation comes 
from (1) the normal-ordering problem in light-front field theory, and (2) the unboundedness 
of the Hamiltonian for the charge conjugation even sector. We think that we can avoid 
these problems if we are confined in the very strong coupling region. We will discuss these 
problems in the final section. 



II. BOSONIZED MASSIVE SCHWINGER MODEL 



A. LFTD for the bosonized massive Schwinger model 



It is well known that the massive Schwinger model [|Tl],|l2| (Q-E-Di+i with massive 
fermion) , 

C = -\f^.F^^' + ^Yi^{id, - eA^) - m]^ , (2.1) 
F^iv = d^Ay — dyA^ , 

has an equivalent bosonic form 

1 

Cb = -df,(pdy - + cm/i : cos(2v^0) : , (2.2) 

where c is numerical constant (c = /2tt, with 7 being the Euler constant), and [i = ej a/vt- 
In the following we set yU = 1, so the strong coupling region corresponds to small fermion 
masses. We quantize this model on the light cone. 

In the equal-time quantization, the normal-ordering is well-defined in the interaction 
picture. (The above normal-ordering is with respect to the free boson with mass yU.) It is 
known, on the other hand, that a naive normal-ordering fails in the light-front quantization, 
because one usually neglects generalized tadpole diagrams [|1^ . The effects of the generalized 
tadpoles amount to the renormalization of the coupling constant . In this paper, however, 
we do not consider this coupling constant renormalization, simply assuming that the effects 
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of it are negligibly small for the very strong coupling region. We will discuss it in the final 
section. 

Expanding (in the Schrodinger picture) in terms of creation and annihilation operators, 



H^-) = f ^^^\-<P^y~''^" + , (2-3) 



[a(p+),at(g+)] - g+) , [a(p+), a(g+)] = [at(p+), at(g+)] = , (2.4) 



we obtain the light-cone Hamiltonian: 



P- ^Ho + V , 

Ho = ^ r%a\p)a{p), (2-5) 
z Jo p"^ 



where /t^ = 1 + Ancm. Note that because we are going to consider the Tamm-Dancoff (TD) 
truncation up to including three-boson states, the interaction terms containing more than 
six creation and/or annihilation operators are irrelevant, and will be ignored hereafter. The 
interaction terms are expressed in terms of and a in the following way, 

cm f'^ dpi + + t 

^4 = — 7r(2^) / n +^^2 +P3 -P4)ai44«4 

4. ^0 Pi 

+6S{pi +P2-P3- P4)ai4«3a4 + 45(pi - P2 - Pa - P4:)a\a2a3a4] , (2.6) 

cm . /■°° dpj J. X ^ J. . 

^ ^ 'rT^ ' / -l-J- — i^^iP^ +^^2 -Fp3 +P4 +P5 -P6)aia2a3«4«5a6 

D. Jo Pi 

-M55(pi P2 + Pa + P4 - P5 - P6)ai444a5a6 
-F205(Pi P2 + Pa - P4 - P5 - P6)aiCi24«4a5a6 
-M55(pi P2 - Pa - P4 - P5 - P6)ai4aaa4a5a6 

-F6(5(pi - P2 - Pa - P4 - P5 - Pe)^! 0203040506] , (2.7) 

where we use the abbreviation = a{pi). 

The spectrum of the model is obtained by solving the Einstein-Schrodinger equation, 
2P+P~ 10) = 10). We truncate the Fock space up to including three-boson states. 
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+ |02)p + |03)p , (2.8) 



= <Pia\V) |0) 

rV 2 



2) p = ^ / n ^'^(Ep^ - V)Hvi,V2)a\p{)a\p,) |0) , 

1 rV dp. 3 

3) p = -j= / n ^'^(Ep« ~ ^)03(pi,P2,P3)a^(pi)aHp2)a^(p3) |0) 

vol -^0 ■_i \/Vi 



where P is the eigenvalue of the momentum operator P+. Note that (f>2{Pi,P2) and 
4>3{pi,P2,P3) are totally symmetric functions of the arguments. 

The boson creation operator can be expressed in terms of the fermion and antifermion 
creation and annihilation operators (see Ref. for the notation), 

CP dh 
a\p)= h\k)d\p-k) 

^0 {2n)^k{p~k) 

POO rl U 

+ / [b^{p + k)b{k)-d\p + k)d{k)]. (2.9) 

^0 {27T)^Jk{p + k) 

It is clear that the symmetry under (p —>■ —(f) (or —a^) is the charge conjugation symme- 
try of the fermionic theory. Because of this symmetry, a state has a definite transformation 
property, 

10); ^10);, 10); -.-10);. (2.10) 

It is easy to see that |0i)p and |03)p are odd while |02)p is even. 

We rescale momenta /V, and the wave functions, 02(pi,P2) and 

03(Pi)P2,P3), are replaced by 02(a;i,X2) and (l)s{xi, X2, x^)/ VV, respectively. The Einstein- 
Schrodinger equation leads to two sets of eigenvalue equations for the wave functions, ac- 
cording to the transformation property under the charge conjugation. The even one involves 
only 02: 



M'^(f)2{xi,X2) = IJ,"^ (— + —) (i)2ixi,X2) - cm{27r) [ ^^i^5(?/i + ?/2 - 
\Xi X2J Jo JmB 



1) 



02(^1,2/2) 



^/mB ^/xiX2 

(2.111 



with xi + a;2 = 1, while the odd one is the coupled equations for 0i and 



2 /■! ^ (ill- 

M20i=/i20i-cm(27r)^ / J] ^^(E " l)'^3(z/i, Z/2, Z/s) , (2.12) 

Vol JO ^^-^ y/yi -^^ 

3 1 

dyidy2r^, , , 03(?/i, ?/2, xi) 



.0 N /"^ dyidy2. , . 
-cm(27r) / + 1/2 + a^i - 1) 



yyiy2 

+ oyVi + y2 + X2- 1) ; h d{yi + y2 + X3~ 1)- 



A 2 f^TT^yijirX^ 03(1/1,1/2, 1/3) 
+ cm(27r)-/ II — (5(2^i/i - 1) 

-nm(27r) ^- , (2.13) 

with xi + X2 + X3 = 1. Note that because the even sector ( |2.11| ) does not depend on Vg, 
the Hamiltonian is not bounded from below. We will see, however, that we get reasonable 
results if we do not employ many basis functions (see below) and are confined in the very 
strong coupling region. 

These complicated equations are converted to two matrix eigenvalue equations by ex- 
panding the wave functions in terms of basis functions. The choice of the basis functions is 
very important for efficient numerical work. We choose the following basis functions so that 
we can calculate the matrix elements analytically. 

MxU^2)=T.Hxi^2y^'^\ (2.14) 

/ 

cj)s{x^,X2,Xs) = Eci(a;iX2X3)'^+'/'(xH4' +4') , (2-15) 
1 

where / = 0, 1, 2, • ■ ■ , A^i, Zi = 0, 1, 2, ■ ■ ■ , A^'2, and /2 = 0, 2, ■ ■ ■ , 2A^3. It is easy to see that any 
symmetric polynomial in xi, X2, X3 with the constraint Xi + X2 + X3 = 1 can be expressed 
by using the above basis functions (up to (xiX2X3)"'^/^) ( p.l5[ ). 



B. States in terms of fermion variables 

How are those states expressed in terms of fermion operators? By using ( p.9|) , it is 
straightforward to express \(p2)-p as follows. 



,(2)\ /"^ dkidk2 

^2 



,(2) 



^ Jo 2'KyJkxk2 



6{k, + k2-V)4'\h,k2)b\4\0), 



1 ^-p 4 ju 4 



^ 2 Jo t=\V^: 



« i=l 



(2.16) 



where 



k2] ks, ki) = 



° V(^i-g)(A:2 + g) 



\/2 



(/Cl + /C3)(/^2 + /C4) 



In a similar way, \4>3)-p is expressed as 



(2.17) 
(2.18) 



dkidk2 



(3) 



5(fci + A;2- (A;i,A;2)6l4|0), 



(3) 



1 /•■P ^ //i- ^ 



-P 2 Jo fj^ 



» i=l 
6 



(2.19) 



^e'O^ = r n -^^m ki - r)4'\h, k2, ks; k,, k,, ke)blblbUldldl |0) , 



■P 3\Jo t=\V^- 



« i=l 



where 



4'\ki,k2)^ 



'ipi\ki,k2;k3,k4 



(/:*/.' 



rki rh (f)^{ki-liM-h,k2 + l2) 

al2- 



V(^i-/i)(/i-/2)(A:2 + /2) 



- dh 
Jo Jo 



Al Mkl-h,k2-l2,ll + l2) , ^ , V 

dh I \ + [ki'^ k2) 



\/{kl-h){k2-l2){h + l2), 



, (2.20) 



^3! r JsHM + h-lM + h) _ ^ _ ^ ^2.21) 



yjl{ki + k^ - I){k2 + k4) 
+ {ki ^ k2, kz ^ k4) 



and 



ilJe\ki,k2,k3;k^,k5,kQ) 



hjh + k4, k2 + fcs, ks + fee) 
(A;i + A;4)(/C2 + A;5)(A;3 + A;6) 



- ik2 ^k3)\ + 



^ ki,k2,ks^ 



\ 



cychc 



(2.22) 



8 



The even state contains the fermion two- and four-body components: 

'r,m)v =V5{V' -V)m + Wl] , 

W^= dxidx25{xi + X2-l)\ijf\xi,X2)\\ (2.23) 
Jo 

„1 4 4 

Wl = ]]_dx,6{J2xi-l)\iljf\xi,X2;x3,x^)\'^ , 



i=l 1=1 

and the odd state contains the fermion two-, four- and six-body components: 



V'{^\4>)v = 


V5{V' - V) [W° + W2 + w°] 






w° = 


0iP + / dXidX2S{Xi + X2 — 

Jo 


l)l^f 


{Xl,X2)\'^ , 


W'l = 


„1 4 4 

/ T\dxiS{y2xi - i)\ipf\xi 

^0 i=l ^=l 




a;4)P , 


^6 = 


»1 6 6 

/ T\dxiS{y2xi - l)\ip^^\xi 

^0 i=l *=1 


X2, 


X4, X5, Xe)^ 



(2.24) 



Once we numerically obtain the eigenstates of the Einstein-Schrodinger equation, i.e., 
the basis function expansion coefficients of the boson wave functions, ( 2.14| ) and ( 2.15|) , we 



can calculate all of these integrals analytically so that we can obtain fermion wave functions. 
This virtue comes from our clever choice of the simple set of basis functions. 

We normalize states in a Lorentz invariant way, ^'(010)^ = V6(V' — V), so that W2 + 
Wl = W2 + W2 + = 1. By saying that a state has a 50.000% fermion two-body 
component, we mean that the state has W2 = 0.50000. 



C. Numerical Results 

1. mass spectrum 

We calculate the invariant masses M for various values of the fermion mass m. We find 
that convergence is good enough for A^^i = 5, A^2 = 4, and N3 = 1 (total number of basis 
functions is 16). The mass spectrum for 0.001 < m < 0.05 is shown in Fig. |^. In the 
following, we concentrate on the case m = 0.01. 
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The lowest state is the meson state. It is charge conjugation odd. It has 99.997% 
one-boson component and 0.003% three-boson component. Its mass is 1.0174. 

The second lowest one should be the two-meson bound state, though its mass appears 
slightly above the threshold, M = 2.0468. It is charge conjugation even. Because of the 
present approximation, its two-boson component is 100%. It may contain a four-boson 
component if we include four-boson states, but we expect that it is negligibly small in the 
strong coupling region. 

The third and fourth states come from the charge conjugation even sector. These are 
regarded as two-meson scattering states. They are completely two-boson states because of 
the present approximation. 

The fifth one is the three-meson bound state, though its mass also appears slightly 
above the threshold, M = 3.0766. It is charge conjugation odd. It has 0.000% one-boson 
component and 100.000% three-boson component. 

The two-meson and three-meson bound states appear below the thresholds when the 
fermion mass gets larger. In the strong coupling region, they should appear just below the 
threshold. (In the strong coupling limit, they are just on the thresholds, not being bound.) 
It is therefore difficult to get them below the thresholds numerically. 



2. Fermion components 

How do these states look like if they are expressed in terms of the fermion variables? 
By using ( p.23|) and ( p.24|) , we get the following results: (1) The meson state is perfectly a 



fermion two-body state (100.00%). (2) The two- meson bound state contains 45.92% fermion 
two-body component and 54.08% fermion four-body component. (Its fermion six-body com- 
ponent is 0% due to the present approximation.) (3) The three-meson bound state, which 
is of our main interest, contains 28.38% fermion two-body component and 63.87% fermion 
four-body component and 7.75% fermion six-body component. 

For the comparison, we give our previous results 0: (1) The meson state has M = 
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1.01813 and is perfectly a fermion two-body state (100.000%). (2) The two-meson bound 
state has M = 2.05612. It has 54.408% fermion two-body component and 45.592% fermion 
four-body component, with httle fermion six-body component. (3) The three- meson bound 
state has M = 3.10814. It has 44.285% fermion two-body component, 53.123% fermion four- 
body component, and 2.592% fermion six-body component. Obviously these are consistent 
with the results obtained by using bosonization, except for the differences in the fermion 
two-body and four-body components of the three-meson bound state. 

The most important result is the demonstration that the three-meson bound state, which 
is almost perfectly (100.000%) a three-boson state, has only 7.75% fermion six-body com- 
ponent. We think that this is a very strong support that the identification we made in the 
previous paper is correct. 

3. Wave functions 

It is interesting to see the wave function of the "relative motion" of the two-meson bound 
state (Fig. It looks very similar to the one we obtained in the previous paper (Fig. 6 
of Ref. [§). We can also calculate the fermion two-body wave function of the two-meson 
bound state 1^2'^ (Fig. ||) which should be compared with Fig. 7 of Ref. 0. The results are 
completely consistent with the previous calculations. 

We next try to figure out how the three-meson bound state looks like by examining the 
wave function. We show |03p of the three-meson bound state in Fig. | (for m = 0.001) and 
Fig. 1^ (for m = 0.05). Because the wave function in momentum space is spread for larger 
fermion masses, and has a sharp peak for smaller fermion masses, we have an intuitive 
picture that for strong couplings the three-meson bound state is loosely bound while for 
weak couplings it has a relatively compact form. We look for asymmetry which indicates 
that two of the three mesons are more closely bound than the third, but we are not able to 
find any. 
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III. SUMMARY AND DISCUSSIONS 



We calculate the mass spectrum of the bosonized massive Schwinger model by using the 
LFTD approximation up to including three-boson states. We showed that the three-meson 
bound state, which is almost 100% a three-boson state, has only a few percent fermion six- 
body component, and the result is consistent with our previous LFTD calculations in terms 
of the fermion variables. We also show that the other quantities are also consistent with the 
previous calculations. By obtaining the wave function of the three-meson bound state we 
are able to have an intuitive picture of the three-meson bound state, namely, it is loosely 
bound for strong couplings while it has a relatively compact shape for small couplings. 

In our present approximation, the two-meson bound state cannot have a non-zero fermion 
six-body component. Of course, if we include four-boson states, it can have a non-zero six- 
body component. We however expect that it will be negligibly small. 

The LFTD approximation for the bosonized theory is good only for strong couplings. 
There are three reasons: (1) The Fock space of the bosonized theory is that of a structureless 
boson, so that we neglect the internal structure of the meson as the first approximation. 
As is known from previous calculations, it is not a good description of the meson to ignore 
the internal structure already at m = 0.1. Note that the nontrivial momentum depen- 
dence of the fermion two-body wave function of the meson, ip2, comes from the many-boson 
components, when expanded in terms of fermion variables. (2) The Hamiltonian for the 
charge conjugation even sector is not bounded from below. (Of course this is because of 
the Tamm-Dancoff truncation.) Because of that we observed instabilities when we increase 
the number of basis functions or when we increase the fermion mass. For small number of 
basis functions and small values of the fermion mass, our method does not scan the 'high 
energy' states and the low energy states are insensitive to the unboundedness of the upside- 
down double well potential. Our choice of basis functions might also suppress the 'decay' 
of the states. (3) The normal-ordering problem becomes serious for weaker couplings. Our 
previous calculation |p shows that the meson mass depends on the fermion mass almost 
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linearly for a wide range of the fermion mass. But the equation (|2.12|) shows, according 



to the variational principle, that the meson mass must be smaller than fi = ^/l + 47rcm, 
which has the linear dependence on the fermion mass only in the strong coupling regions. 
For weaker fermion masses, this appears to put a stringent "upper bound." Of course this 
should not be true. (The structureless one-boson approximation to the meson cannot be 
better than that including up to fermion six-body states, especially for weak couplings.) The 
"coupling" c must be renormalized and must have a nontrivial fermion mass dependence. 



The investigation in this direction is now in progress ||T^. For the purpose of the present 
paper, however, it is sufficient to notice that the renormalization effects are negligibly small 
for small fermion masses. The fermion mass dependence is almost identical to that of our 
previous calculation. 

Another pathology of the bosonized model can be seen in the wave function of the 
meson state. It is known that the fermion two-body wave function of the meson state, 

',(3) 



'ip2{x, 1 — x) = 01 + ^/'2 {x,l — x), must vanish at x = and x = 1 in the massive theory [jl6 
The calculated wave function shown in Fig. ^ does not satisfy this requirement. Furthermore, 
it is not a concave function of x. (See Fig. 5 of Ref. [0, for example.) We do not know why 
it does not have the correct behavior. 

It would be interesting to use the bosonized theory to investigate the effects of the theta 
vacuum 



1^, though the investigation in the original fermionic theory requires an intensive 



study O . 
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FIG. L The mass spectrum for 0.001 <m< 0.05, at Ni = 5,N2= 4, and A^s = L The dashed 
and dotted hne stand for the two-meson and three-meson thresholds respectively. 
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FIG. 2. Squared wave function of the relative motion of the two-meson bound state, 
\(p2ix, 1 — is shown for various values of the fermion mass. This should be compared with 



Fig. 6 of our previous work. 




X 



FIG. 3. Fermion two-body wave function for the two-meson bound state, ^p2 is shown 

(solid line) with that obtained by the previous calculation in terms of fermion variables (dashed 
line) for m = 0.01. 
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FIG. 4. Squared wave function of the three-meson bound state, \(j)3{xi,X2,X3)\'^, for m = 0.001 
projected on the xi-X2 plane. Only the region (xi > 0, X2 > 0, xi + X2 < 1) is the support. 




FIG. 5. Squared wave function of the three-meson bound state, |(/)3(xi, X2, xa)^, for m = 0.05 
projected on the X1-X2 plane. Only the region (xi > 0, X2 > 0, xi -|- X2 < 1) is the support. 
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FIG. 6. Fermion two-body wave function of the meson, \ip2{x,l — x)p is shown for various 

values of the fermion mass, where iIj2 = (pi + tp^'^ ■ This should be compared with Fig. 5 of our 

previous work. 



19 



